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Abstract 

We propose a scale-free network model with a tunable power-law 
exponent. The Poisson growth model, as we call it, is an offshoot 
of the celebrated model of Barabasi and Albert where a network is 
generated iteratively from a small seed network; at each step a node is 
added together with a number of incident edges preferentially attached 
to nodes already in the network. A key feature of our model is that 
the number of edges added at each step is a random variable with 
Poisson distribution, and, unlike the Barabasi- Albert model where this 
quantity is fixed, it can generate any network. Our model is motivated 
by an application in Bayesian inference implemented as Markov chain 
Monte Carlo to estimate a network; for this purpose, we also give a 
formula for the probability of a network under our model. 



Keywords: Bayesian inference, complex networks, network models, power- 
law, scale-free 

1 Introduction 

Until recent times, modeling of large-scale, real- world networks was primar- 
ily limited in scope to the theory of random networks made popular by |Erdos 



and Renyi (1959). In the Erdos-Renyi model, for instance, a network of N 



nodes is generated by connecting each pair of nodes with a specified probabil- 
ity. The degree distribution p{k) of a large-scale random network is described 
by a binomial distribution, where the degree A; of a node denotes the number 
of undirected edges incident upon it. Thus, degree in a random network has 
a strong central tendency and is subject to exponential decay so that the 
average degree of a network is representative of the degree of a typical node. 

Over the past decade, though, numerous empirical studies of complex 
networks, as they are known, have established that in many such systems — 
networks arising from real-world phenomena as diverse in origin as man- 
made networks like the World Wide Web, to naturally occurring ones like 
protein-protein interaction networks, to citation networks in the scientific 



Albert et al. 


(1999 


), 


Jeong et al. 


(2001 


), and Redner 


(1998 



respectively — the majority of nodes have only a few edges, while some nodes, 
often called hubs, are highly connected. This characteristic cannot be ex- 
plained by the theory of random networks. Instead, many complex networks 
exhibit a degree distribution that closely follows a power-law p{k) oc k~'^ 
over a large range of k, with an exponent 7 typically between 2 and 3. 
A network that is described by a power-law is called scale-free, and this 
property is thought to be fundamental to the organization of many complex 
systems; Strogatz (2001). 



As preliminary experimental evidence mounted (see Watts and Strogatz 



(1998), for instance), a simple, theoretical explanation accounting for the 



universality of power-laws soon followed; the network model of Barabasi and 



Albert ( 1999 ) (BA) provided a fundamental understanding of the develop- 



ment of a wide variety of complex networks on an abstract level. Beginning 
with a connected seed network of to nodes, the BA algorithm generates a 
network using two basic mechanisms: growth, where over a series of itera- 
tions t = to,to -\- 1, . . . the network is augmented by a single node together 
with m < to undirected incident edges; and preferential attachment where 
the new edges are connected with exactly m nodes already in the network 
such that the probability a node of degree k gets an edge is proportional 
to r{k) = k, the degree of the node. When m is fixed throughout, Bol- 



lobas et al. (2001 ) showed rigorously that the BA model follows a power-law 
with exponent 7 = 3 in the limit of large t. The sudden appearance of 
the BA model in the literature, nearly a decade ago, sparked a flurry of 
research in the field, and, consequently, numerous variations and generaliza- 
tions upon this prototypal model have been proposed; [Albert and Barabasi 



(2002) and Newman (2003) rank as the preeminent survey papers on the 
subject. 

In this paper, we propose a new growing model based on preferential 
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attachment: the Poisson growth (PG) modeL Our model, as described in 
Section |2] is an extension of the BA model in two regards. Firstly, we 
consider the number of edges added at a step to be a random quantity; 
at each step, we assign a value to m according to a Poisson distribution 
with expectation A > 0. Secondly, we avail ourselves of a more general 
class of preferential attachment functions r{k) studied by several authors 
including [Krapivsky and Redner (20011 and Dorogovtsev et al. (2000). In 
Section |3] we argue that the degree distribution of the PG model follows a 
power-law with exponent 7 that can be tuned to any value greater than 2; the 
technical details of our argument are left for the Appendix. In addition, we 
conducted a simulation study to support our theoretical claims. Our results, 
provided in Section |4j show that the values of 7 we estimated from networks 
generate under the PG model are in agreement with those predicted by our 
formulae for the power-law exponent. 

Our motivation for proposing the PG model, as explained in Section [5] 
arises from a need for a simple, yet realistic model that is serviceable in ap- 
plications. In fact, with our model every network has a nonzero probability 
of being generated, in addition to possessing a tunable power-law exponent. 
In contrast, the BA model has a fixed 7, and is subject to numerous struc- 
tural constraints which severely limit the variety of generable networks. We 
give a simple formula for the probability of a network under the PG model, 
which can be applied quite naturally in Bayesian inference using Markov 
chain Monte Carlo (MCMC) methods. Firstly, given a network G, we may 
estimate the PG model parameters, or engage in model selection in the case 
when we have more models; or, going against the grain, we may estimate 
an unknown G from data using our PG model formula as a scale- free prior 
distribution. 

Finally, other scale-free models have been put forth in the literature, 
besides our new model, that are realistic enough for use in applications. An 
extension of BA model by Albert and Barabasi (20001 incorporates "local 
events," which allows for modifications, such as rewiring of existing edges, 
to the network at each step. Other authors, including iSole et al. (20021, pro- 
posed a class of growing models based on node duplication and edge rewiring. 
Other scale-free model not based on growth have also been proposed; for ex- 
ample, the static model of Lee et al. (20051. Among these models, the PG 
model is the simplest preferential attachment model sufficiently realistic for 
use in applications. 
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2 The Poisson Growth Model 



In the PG model, we begin with a small seed network of to nodes. Let 
Gt = {Vt,Et) be the network at the onset of time step t > to where Vj = 
{vi,V2, ■ ■ ■ ,vt} is a set of t nodes and Et is a multiset of undirected edges 
so that multiple edges between nodes, but no loops, are permitted. The 
updated network Gt+i is generated from Gt as follows: 

Poisson growth: A new node ff+i is added to the network together with 
rrit incident edges; nit is a random variable assigned according to a 
Poisson distribution with expectation A > 0. 

Preferential attachment: Each edge emanating from vt+i is connected with 
a node already in the network. Node selection can be considered as a 
series of mt independent trials, where at each trial the probability of 
selecting a node from Vt with degree k is 

Qtik) = —-^ , (1) 

where ki^t is the degree of node Vi at step t. Define Si^t as the number of 
times node Vi is chosen at step t. Then the entire selection procedure 
is equivalent to drawing a vector {si^t, S2,t, ■ ■ ■ , st^t) from a multinomial 
distribution with probabilities qt{ki^i), . . . , qt{kt,t) and sample size mt. 
Equivalently, Si^t has Poisson distribution with expectation Xqt{ki^t) 
independently for i = 1, . . . ,t. 

The PG model is determined by the choice of r{k); we concentrate on 
two specifications and discuss their implications in the next section. Firstly, 
let 

r{k) = k + a (2) 
where the offset a > is a constant. More generally, wc define 

r{k) = k + a, k>l, and r(0) = b (3) 

by taking o > — 1 with extended domain, but in doing so define a threshold 
parameter 6 > 0. Indeed, the latter formulation includes the former as a 
special case when we restrict a = 6 > 0, so that overall our model is specified 
by the parameter 6 = (a, b, A) . 

The BA model can be explained as a reduction of our model by taking 
= 6 = 0, and by fixing 1 < = m < to so that the number of edges added 
to the system at each step is a constant; the new edges are preferentially 
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attached from the new node to exactly m other nodes. Many structural 
constrains are imphcit in the BA modeL Indeed, at step t, a network with t 
nodes must have m{t — to) + \EtQ \ edges, none of which are multiple, whereas 
the number of edges for the PG model can take other values. 

A number of extensions of the BA model based upon generalizing r{k) 



have been proposed. In particular, Krapivsky et al. (2000) analyzed a version 



where the preferential attachment function is not linear in the degree A; of a 
node, but instead can be a power of the degree k'^, i' > 0. They showed that 
for the scale-free property to hold, r(k) must be asymptotically linear in k. 



In a subsequent work, Krapivsky and Redner (2001) and Dorogovtsev et al 



(20001 independently went on to establish that adding the offset a > —m as 
in ([2j) does not violate this requirement, and derived the power-law exponent 
7 = 3+a/m. Their result is analogous with our reported power-law exponent 



in (|7[) w ith A = m as seen in the next section. Furthermore, Krapivsky and 
Redner (2001 ) investigated an attachment function similar to ([3]) defined by 



r{k) = k,k > 2,r(l) = 6, r(0) = 0. As they took m > 1 they did not need 
to be concerned with nodes of degree k = 0. The power-law exponent they 
derived in this case is reminiscent of our result in ^ . 

3 The Degree Distribution of the Poisson Growth 
Model 

In this section, we discuss the degree distribution p(k) for networks generated 
under the PG model. The main result is that the degree distribution follows 
a power-law 

p{k) ~ k-\ (4) 

where ~ bk indicates these two sequences are proportional to each other 
so that ak/bk converges to a nonzero constant as /c — > oo. This result is an 
immediate consequence of the recursive formula 

{k + a-l+-f)p{k) = {k + a-l)p{k-l) (5) 

for sufficiently large k, and thus p{k) {k + a — 1)^^ k^'^ . The power-law 
exponent is 

^^^^a + ib-a)pm 
A 

for the preferential attachment function defined in (|3]), and the exponent 
takes the range 7 > 2; the lower limit 7 — > 2 can be attained by letting 
a = — 1, 6 = 0, and A — > 0. This lower limit is in fact the limit for any form 
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of r{k) when A does not depend on t; 7 must be larger than 2 to ensure 
the mean degree ^^=Qkp{k) = 2A converges. For the special case ([2]), the 
exponent becomes 

7 = 3 + p (7) 

and the range is 7 > 3. 

To make the argument precise, we have to note that Gt is generated 
randomly and the degree distribution of Gt also varies. Let nt{k) be the 
number of nodes in Gt with degree k. Since X]fc^o^*(^) ~ observed 
degree distribution of Gt is defined by pt{k) = nt{k)/t for A; > 0. For 
sufficiently large t, pt{k) may follow the power-law of Q as seen below. 

We consider a moderately large k for the asymptotic argument as t ^ 00. 
The maximum value of k for consideration \s k ^ for a given t with a 
constant c = 1/(7 + 2 + e) with any e > 0. Then, the expectation of pt{k) 
can be expressed as 

^ptm ~ k-\ (8) 

which is the power- law we would like to show for pt{k). The variance of 
Pt{k) will be shown as 

\{pt{k)) = 0{e+H-^), (9) 

indicating the variance reduces by the factor 1/t. Note that ^ is not a tight 
upper bound, and the variance can be much smaller. See the Appendix for 
the proof of (|8| and ([o]). Let < d < 1/(27-^2 -he), and consider k = 0(t''), 
which is even smaller than t'^. Then, 

"^^^jf^^ = 0[k^^^^^'H-^l^) = 0{n (10) 

with a = ci(7 -|-l-|-e/2) — 1/2 < 0, and thus the limiting distribution 
limt^oo Pt{k) = p{k) follows the power-law of Q. By taking e — > 0, the 
power-law of pt{k) is shown up to A; ~ with d < 1/(27 + 2). 

It remains to give an expression for p(0) in ([6|. We will show in the 
Appendix that p(0) is a solution of the quadratic equation 

(6 - a)x^ + {2\ + a + \h-{h- a)e-^)x - (2A + a)e'^ = 0. (11) 

For a^h^ one of the solutions 

p(0) = ^ ^_ [{(2A + a + Ab - (6 - a)e-^f + 4(6 - a)(2A + a)e-^}^^^ 



-(2A -h a A6 - (6 - a)e 



fl2) 
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is the unique stable solution with < p(0) < 1; this can be checked by 
looking at the sign of pt_|_i(0) — pt{0) in the neighborhood of p(0). 



4 Simulation Study 



A small simulation study was conducted to support our theoretical claims 
of Section |3j Specifically, we wish to confirm via simulation that the degree 
distribution p{k) of Q as well as its expected value E(pt(A;)) as in ^ follow 
a power-law with 7 as in ([6]). To that end we generated networks under the 
PG model for a variety of parameter settings. For each specification of 9 we 
generated n^jm = 10^ networks of size N = 5000, each from a seed network 
of a pair of connected nodes. We included the BA model, generated under 
analogous conditions, so as to demonstrate the soundness of our results 
which are summarized in Table [H 

In point of fact, estimating 7 from a network can be quite tricky and 
it has been the subject of some attention in the literature; see |Goldstem 
et al. (2004). We sided with the maximum likelihood (ML) approach de- 



scribed by Newman ( 2005 1 . In this methodology, the ML estimate of 7 for 
a particular network is given by 



-1 



7 = 1+ Yl \ \Y. ^(^)iog 



k 



where n{k) is the number of nodes with degree k, and kmin is the minimum 



degree after which the power-law behavior holds. Bauke (20071 studied 
selecting a value for kmin by using a goodness of fit test over a range of 
kmin] however, we shied away from this level of scrutiny as we found that 
taking kmin = 10 was reasonable for our examples. This methodology is 
illustrated in Figure [T] (a) and (b) where we plot the degree distribution with 
7 for a typical network generated by the BA and PG model, respectively. 

Returning to Table [T in each case, we confirm (j4]) has power-law ex- 
ponent as predicted by (6| and (12). We computed 7 for each network, 
and calculated the mean and standard deviation of 7 values for Usim net- 
works. We observe that the mean 7 agrees well with the predicted 7, and 



the variation of 7 is relatively small as suggested by (10) 



In addition, to show that the same holds for ([8|), in each case we com- 
puted the average degree distribution of the Usim networks as an estimate 
of E(p(A;)). Then we estimated the degree exponent ^avg as seen in the table 
and Figure |2j Again, the simulated results match well with theory. 
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5 Discussion 



The PG model has a special place in the class of preferential attachment 
models. It has a tunable power-law exponent and a simple implementation, 
yet it can generate any network. In contrast, the BA model and its gen- 
eralizations described in Section [2] have serious restrictions on the types of 
networks that can be generated because m is held constant. For example, 
at step t an instantiation of the BA model will consist of a t node network 
with the number of edges equal to exactly m(t — to), plus the number of 
edges in the seed network. The simple design of our model makes comput- 
ing the probability of a network straightforward. This in combination with 
its modeling potential gives rise to several useful applications in Bayesian 
inference. 

In explicit terms, let G = (V,E) be a network with = nodes where 
V = {v'l, . . . , v'j^}. Furthermore, let Gat = (Vat, En) be a network generated 
under PG model after step A — 1 so that Vjy = {vi, . . . , f at}, where the seed 
network consists of a single node. The association between V and Vjsr is 
defined by a permutation a = {ai, . . . , ajy) so that Vi = v'^_. Given G, once 
we specify a, then it is straightforward to compute ki^t, Si^t for i = 1, . . . ,t; 
t = 1, . . . ,N — \. Then the probability of G given 9 = (a, b, A) and a is 

P(G\9,a) = Yl |^ne"^'^'(^-)^^^fe^^ . 

One application is when G is known and we wish to estimate 9. This 
can be done by assigning a prior 7r(^) for and the uniform prior on a. The 
posterior probability of {9, a) given G is 

7T{9,a\G) ocF{G\9,a)n{9). 

Using MCMC to produce a chain of values for {9, a), the posterior 7r(6'|G) 
is simply obtained from the histogram of 9 in the chain. Moreover, this 
procedure can be used for model comparison, if we have several models for 
generating the network. 

Another application is when we wish to make inference about G from 
data D with likelihood function P{D\G). The posterior probability of 
(G, 9, a) given D is 

TT{G,9,a\D) oc P{D\G)F{G\9,a)7r{9). 

Then the posterior 7r(G|-D) is simply obtained from the frequency of G in 
the chain. Indeed, we used this approach for inferring a gene network from 



microarray data in Sheridan et al. (2007). 
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Recall that the PG model produces networks with multiple edges. In 
practice, we often want to restrict our interest to networks without multiple 
edges. As an approximation, we could apply the formula for P{G\9, a) just as 
well in this case. Alternatively, we propose a slight modification to our model 
where we generate nit edges at step t according to a binomial distribution 
with parameter p = X/t and sample size t. In this formulation the seed 
network must be selected such that X < to, otherwise p > 1 may occur. 
Then by sampling nodes without replacement, multiple edges are avoided. 
In our simulation (results not included) we found that these modifications 
do not change the power-law. 

Finally, though we made specific choices for r{k) in our arguments, the 
PG model can be generalized to a wider class of preferential attachment 
functions. For instance, Dorogovtsev and Mendes (2001) investigated accel- 
erated growth models where mt increases as the network grows. It should 
be possible to incorporate accelerated growth into PG model by gradually 
increasing the value of A over time. Another line of generalizations of the 
PG model is via the inclusion of local events. 



Appendix: Proofs 



The expected value of pt{k) 

Here we give the proof of ([s]). We assume that the functional form of 
r{k) is and a modification to handle ([s]) is mentioned at the bottom. 

Let I{A) denote the indicator function of the event A, so I{A) = 1 A 
is true and I{A) = if ^ is false. We use the notation P(-), E(-) and V(-) 
for the probability, expectation and the variance, and also P(-|^), E(-|74) 
and V(-|^) for those given a condition A. By noting 

t 

nt+i{k) = ^ I{ki^t + Si^t = k) + I{mt = k), 

i=l 

the conditional expectation of nt+i{k) given Gt is 

t 

E{nt+i{k)\Gt) = Y,P{ht + Si,t = k\Gt)+F{mt = k\Gt) 

i=l 
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The last term e~^X^/kl ~ {e\/k)^ can be ignored for a large since it 
is exponentially small as k grows. We examine the terms in the summation 
over s = 0, 1, . . . , A; for k = 0{t^) as t ^ oo. For a fixed s, qt{k — s) ~ k/t for 
a linear preferential attachment model. More specifically, for r{k) = k + a, 
k>0, 

/, N rik — s) k — s + a. ^, i/o,, 
qt{k -s)= ^ ' = (1 + 0{t-^/^)), 

because the mean degree of Gt is 



t-i 



t'=to 



and the denominator of qt{k) is 
t t 

r{ki,t) = Y.{ki,t + a)= t{2X + a + 0{t-'/^)). (14) 



i=l 1=1 



Thus the sum in (13) over s = 0, 1 becomes 

M«.- + <i) , ^. [ Mk + a-r, „,, ,-3/2. 



For s > 2, each term is ~ nt{k — s){k/t)^. By noting Yl'^=2^t{k — s) < t, 
the sum over s = 2, . . . ,k becomes 0{k'^t~^). 

Next, we take the expectation of (13) with respect to Gt to obtain the 
unconditional expectation E(nt+i{k)), and replace nt{k) = tpt{k). Using 
the results of the previous paragraph, we get 

nPt+i{k)) = nPt{k)) - ^^XT^ ((^' + ^ + 0{kt-'l^))nPt{k)) 

-{k' + 0{kr^'^))E{pt{k - 1)) + 0{kh~^)\ (15) 



with k' = k + a — 1 and the 7 of ([t]). Let us assume E(pt(/i; — 1)) ~ (k — 1) 
and remember c < 1/(7 + 2). By taking the limit t ^ 00 and equating 
E{pt+i{k)) = E{pt{k)), we get 

(A:' + 7 + o{l))E{pt{k)) = ik' + o{l))E{pt{k - 1)). 
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So that, for sufficiently large t, 

E{pt{k)) ~ k"' 

also holds for k. Since E(pt(A;)) = 0(1) for a fixed k, the power-law holds 
for any k by induction up to k ^ f^. 

For r{k) of the preferential attachment is modified to 

r{k) = k + a + {h- a)I{k = {)), k>0. 



This changes the the denominator of qt{k) in (14) to 
t 



rih^t) = t(2X + a + {b- a)pt{0) + 0{t-'/^)^ , (16) 



i=l 



and thus 2A + a in the updating formula ( 15 ) is replaced with 2A + a + (6 ■ 



a)p(O), leading to Note that pt{0) = p{0) + 0{t~^/^) from ^ shown in 
the next section. 

The variance of pt(k) 

Here we give the proof of ([9| by working on Y{nt{k)) = t^\{pt{k)). 



Although r{k) of ^ is again assumed, the argument is basically the same 
for ([3]). By noting the identity 

\{nt^i{k)) = E(V(nt+i(fc)|Gt)) + V(E(nt+i(fc)|Gt)), (17) 

we evaluate the two terms on the right hand side. 

The conditional variance of nf+i(A;) given Gt is evaluated rather similarly 
as the conditional expectation of (13). By noting \{I{A)) = P(^) — P(A)^, 
\{nt+i{k)\Gt) is expressed for k = 0{t'^) as 

^ ^ ^-\qt(k-s) {Mt{k - S)Y _ ( Xqt(k-s) {Mt{k - S)Y ^ ^ 



s=0 



s\ \ s\ 



,,.\{k + a) ,A(A;-|-a — 1) , , 

-nAk)--\ f + nt(A;-lH^- (18 

^ '{2\ + a)t ^ ' {2\ + a)t ' ^ ' 

where terms from I{mt = k) are ignored for a large k. Thus, the first term 
in (17) is 

E(V(nt+i(A;)|Gt)) =0(A;-^+i). 
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On the other hand, the second term in (17) is evaluated by considering 
the variance of ( 13 ) as 

\{E{nt+i{k)\Gt)) < V(n,(fc)) f 1 - + + 0(fcr^/2) 



{2X + a)t 

+2VnMWwnnt{k - 1)) (^1^^ + o{kt-y')^ 

+V{nt{k - l))0{kH-^) + y/V{nt{k))0{kh-^) + ©(fcV^). 



We substitute these two expressions for those in (17). We will show, by 
induction, that 

Y{nt{k)) < Ak^+H (19) 

holds for all (t, k) with k = O(t^) using some constant A. Let us assume 
that (19) holds for (t, k) and {t, k — 1). By taking a sufficiently large A, we 
have 

V(nt+i(A:)) < Ak'^+'{t - (2A + a)"^) + o{k^+'^^) < Ak'^+'{t + 1), (20) 



implying that (19) also holds for {t + l,k). 

On the other hand, for any random variable < n < t with its ex- 
pectation E(n) fixed, the largest possible variance 0(t)E(n) is attained if 
the probability concentrates on the extreme values and t. Applying this 
upper bound to nt{k) with k ~ f^, we obtain \{nt{k))/t = 0{E{nt{k))) = 
0{k-^t) = 0(A;2+^), implying that ^ holds for any (t, k) with A: ~ t". 

For induction with respect to k, we only have to show 



\{nt{k)) < v{k)t 



(21) 



for a sufficiently large k so that terms from I{mt = /c) in ( 18 ) can be ignored. 
v{k) is an arbitrary constant depending on k. We start from = 0. First 
note that 

t 

nt+i(0) = ^iht = n s^,t = 0) + I{mt = 0). 

i=l 



Thus E{nt+i{0)\Gt) = nt(0)e-^''*(°) + e"^, and so 



V(E(ntH 



Gi))=V(nt(0)) 1 



2Aa 



(2A + a)t 

On the other hand, V(nt_|_i(0)|Gt) is expressed as 



+ 0{t 



-3/2n 



nt(0)(e-^«'(0)-e-2^'''(°))+, 



+ 2nt(0)(l -e~^'?*(°))e' 
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By substituting these two expressions for those in (17), we observe that the 



increase of the variance, i.e., V(ni+i(0)) — V(n((0)) is bounded by a constant, 
and we have V(nt(0)) = 0{t). 

Let us assume (21 ) holds up to A; — 1. Then V{nt+i{k)) can be expressed 



quite similarly as (20), but E(V(nt+i(fe)|Gt)) includes additional terms from 



limt = k); \{I{mt = k)) = 0(1) and E(^^^, CoY{I{ki,t + Si,t = k),I{mt = 
k)\Gt))- For ki^t = ^, the covariance term < P(mt = /c)(l— P(sj,f = 0|Gt)) = 
0(t~^), and for ki^t = k — s with s > 1, the covariance term < P(sj^t = 
s)(l — P(mt = k)) = 0{t~'^). Thus, by taking the sum over i = 1, . . . , it 
becomes 0{t-t'^) = 0(1). Therefore, V{nt-\-i{k)) — V{nt{k)) is bounded by 



a constant, and (21) holds for k. By induction, (21) holds for any k. 



The equation of p{0) 

Here we derive ( |11[ ) for the r{x) of ([3|. By taking the expectation of 
E{nt+i{0)\Gt) = nt(0)e-^^'(°) + e'^ with respect to Gt, and using ([l6|, we 
get 



E(nt+i(0)) =E(nt(0)) 1 



Xb 



{2X + a + {b - a)p{0))t 



+ 0(t-3/2) ) +e 



By substituting nt{0) = tpt{0) and taking the limit t — > oo, we get a formula 
for f{x) = {t + l)(E(pt+i(0)) - E{pt{0))) as a function of x = p{0) 



x{l + 



Xb 



2A + a + (6 - a)x 



+ e 



The quadratic equation (11) is obtained by letting /(x) = 0. In addition, 
the condition df{x)/dx < was checked for the stable solution. 
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Table 1: Summary of estimated power-law exponents from simulated net- 
works. The last column is theoretically predicted 7. 



Model 


Parameters 


Mean k 


Mean 7 it s.d. 


^avg 


7 


BA 


m = 1 


2.0 


3.03 ±0.15 


3.03 


3 


PG 


e = (0,0,1) 


2.0 


3.03 ±0.12 


3.03 


3 


PG 


^= (-0.9,0.1,1) 


2.0 


2.54 ±0.10 


2.51 


2.44 


PG 


= (-0.9,0.1,3) 


6.0 


2.86 ±0.05 


2.86 


2.72 


PG 


e = (0.5,0.5,3) 


6.0 


3.15 ±0.05 


3.15 


3.17 
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Figure 1: Degree distribution p(k) of a typical network plotted on a log-log 
scale with the power-law line using estimated exponent 7. (a) Generated 
under the BA model; 7 = 3.03. (b) Generated under the PG model with 
e = (0,0,1); 7 = 3.01. 
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Figure 2: Avrage degree distribution E(p(fc)) of the simulation with the 
power-law line using estimated exponent "favg- Ploted for (a) the BA model 
and for (b) the PG model with 6 = (0,0,1), where ^avg = 3.03 for both 
cases. 
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